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PROBLEMS IN THE THEORY OF ORDINARY LINEAR DIFFER- 
ENTIAL EQUATIONS WITH AUXILIARY CONDITIONS 
AT MORE THAN TWO POINTS* 

BY 

CHARLES E. WILDER 

In this paper is considered a differential system consisting of an ordinary 
linear differential equation and auxiliary conditions involving linearly the 
values of the solution and its derivatives at interior points, as well as at the 
end points, of the interval over which the equation is considered. The notions 
of Green's function and adjoint systems, already introduced by Birkhofff 
and further developed by B6cherJ for auxiliary conditions involving the end 
points only, are extended to the more general system. The problem of the 
expansion of arbitrary functions in terms of the characteristic solutions of 
the system is stated, but the convergence proof is given in another paper. § 

The author wishes to express his gratitude to Professor B6cher at whose 
suggestion he undertook the work, and to Professor Jackson under whose 
guidance most of it has been done. 

1. Thk auxiliary condition pboblem 
Consider the differential expression 

and the equation L{u) = p{x) , in which p{x) , pi(x) , •••, pn{x) , are 
functions of the real variable x , continuous together with their derivatives of 
all orders in the closed interval ( o , b) . Let ai = a, a^, as, • • • , Ofc = 6 , be 
k points of this interval arranged in order of ascending algebraical magnitude. 
If ^ ( a; ) is any function of x which has at these points derivatives of the first 

* Presented to the Society, April 29, 1916, as part of a more extensive paper, with a slightly 
different title; the remainder of the paper has already been published in these Trans- 
actions, vol. 18 (1917), pp. 415-442. 

t Birkhoff, these Transactions, vol. 9 (1908), pp. 373-395. 

t Bdcher, these Transactions, vol. 14 (1913), pp. 403-420. 

§ See first footnote above. 
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n — 1 orders, then let 

Wij (4,) = aijct>iaj) + a'ij 0' (ay) + oc'/j ct>" {aj) + ■ ■ • + a^tf'^ (A^""" (ay) 

ij = 1,2, ■■•,k), 
and let 

in which an, a'lj , a'/j , • • • , a%~^\ are constants, real or complex. 
If IT,- represents a real or complex number, then the system 

L{u)=p, Wi(u) = Ti (i = 1,2,3, ••-,«), 

consisting of the differential equation with n auxiliary conditions, is called the 
complete system; the system 

(1) L{u)=p, Wi(u)=0 (i = 1,2,3, ••-,«), 
is called the semi-homogeneous system; and finally the system 

(2) L{u)=0, Wi(u)=0 (i = 1,2,3, •••,n), 

is called the reduced system. 

The reduced system is said to be compatible if it has solutions other than 
zero, otherwise it is incompatible. It has r-fold compatibility if it has r and 
only r Unearly independent solutions.* By direct substitution of the general 
solution of L{u) = in the auxiUary conditions Wi{u) = it is found that 

IfyiiVi! ' ' ' ,yn,is a fundamental system of solutions of L{u) =0,a neces- 
sary and sufficient condition that the reduced system haw r-fold compatibility is 
that the following determinant be of rank n — r: 

(3) D ^\Wiiyj)\ (i,y = l,2, ...,n). 

If we denote by uq a particular solution oi L (u) = p , then it is found in 
the same way that 

If the reduced system has r-fold compatibility, a necessary and sufficient con- 
dition that the complete system have a solution is that the matrix obtained by 
bordering the determinant D on the right by the quantities Wi (uo) — iTi, i = 1 , 
2 , • • • , n, be of rank n — r . 

Whence, in particular, for r = , 

If the reduced system is incompatible, the complete system always has a solution 

That this solution is unique follows at once from the fact that the difference 
of any two solutions of the complete system is a solution of the reduced system. 

Assuming that the auxiliary conditions are linearly independent, we shall 
now define a function, Gix, s) , called the Green's function for the reduced 
system, by requiring that for every value of s in the interval ( a , b) , it shall 
have as a function of x the following properties: 

* B6cher, loc. cit., p. 405. 
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1." It is continuous in the interval (a, b) , together with its first n — 2 
derivatives; 

II. It satisfies Xr(M) = at all points of the interval except the point s, 
where its ( n — 1 )th derivative has a finite jump of magnitude 1 ; 

III. It satisfies the auxiliary conditions WiiO) = 0, i = 1, 2, • • • , n. 
The proof that these properties completely define the function is as follows : 
Since G(x, s) satisfies L(u) = , it can be written in terms of a funda- 
mental system of solutions, 

Gix,s) = ci(s)yi(x) + C2(s)y2ix) + ••• + c„ («) j^„ (a;) , x ^ s; 

= di{s)yi(x) + d2(s)y2ix) + •■• +dn(s)yn{x), x^s. 

Now 6(x, s) and its first n — 2 derivatives are continuous at the point s , 
while the ( ra — 1 )th derivative has a finite jump of magnitude 1 there, whence 

ci(«)2/i(«) + ••• + c„(5)2/„(5) = ix{s)yx{s) -\- ■•■ + dn(s)yn{s), 

ciis)y'iis) + ••■ +Cnis)y'nis) = di(s)y'i{s) + ■•■ + dnis)y'„is) , 

Ciis)yr'Hs) + ••• +cnis)yr'\s) = diis)i/r''is) + ••• 

+ dnis)y':-'\s), 
ci{s)y^r'\'i) + ■■■ +c„(*)2/r"(*) = dUs)j/r'Hs) + ••• 

+ d«(5)2/<r-"(«) -1- 

If we set di (5) — Ci(s) = Zi(s) , we can write these equations in the form 

n 

Jlziis)y^fis) = (y = 0,1,2, ...,n-2), 

1=1 

t,Ziis)y<r"{s) = 1. 

The determinant of these equations is the wronskian of the y's and so is not 
zero, since they are linearly independent. Hence the z's are uniquely deter- 
mined and are continuous functions of s in the interval (a, b) . Knowing 
the 2's, we determine the c's and the d's by means of III. If we assume for 
purposes of proof that s is in tlie interval (oi, a^) , we have 



WiiG) ^i:(cjj: Wi,{yi) + dj t, WiAyj)) = 

^■=1 \ <;=1 q=m J 

( i = , 1 , 2 
and on substitution of dj — Zj for c, this becomes 
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These are n equations in the n unknowns, the d's, and their determinant iS D , 
so that if D does not vanish the unknowns are uniquely determined, and then 
by the relation d, — 2!y = cy the c's are also. Hence if the reduced system is 
incompatible it possesses one and only one Green's function. 
By direct substitution of 



U{x}= f pis)Gix,s)> 

Ja 



in the system (1) it is found that the function U {x) is the solution of the 
semi-homogeneous system, by reason of the three defining conditions for 
G{x,s).* 

If Pi has derivatives of the first n — i orders then the first n derivatives 
of the functions Zi (called adjoint functions) exist, and satisfy the differential 
equation 

called the adjoint equation.^ From this, and the explicit expressions for 
G{x, s) in terms of the 2/'s, we can deduce easily the following properties of 
G{x, s) as a function of s : 

I. It is continuous in the interval ( a , b) , together with its first n — 2 
derivatives, except at the points a,- , where it has finite jumps whose magnitudes 
are in general functions of x; 

II. It satisfies M {v) = a,t all points of the interval ( o , b) , except the 
above, and the point x, where its (n — l)th derivative has a finite jump of 
magnitude ( — 1 )"'"^ . 

Since G {x, s) as a, function of s has discontinuities it is apparent that an 
adjoint system in the accepted sense of the termj does not exist. In con- 
sidering the expansion problem it becomes necessary to have functions analo- 
gous to the solutions of the adjoint system in the two point case. These 
may be obtained as the solutions of the adjoint to the integral equation which 
is equivalent to the differential system, but it seems worth while to indicate 
briefly how they may be found without the aid of the theory of integral equa- 
tions. 

The fact that the Green's function as a function of its second argument 
has discontinuities at the points at suggests that it might be profitable to 
investigate the case where it has similar discontinuities with regard to its 
first argument. Such a Green's function would be obtained if the auxiliary 

* It is to be noticed that down to this point the operators Wi might have been any linear 
operators, so that the analysis has a considerably wider range of applicability than has been 
indicated. 

t Schlesinger, Handbuch der Theorie der linearen Differentialgleichungen, vol. I, p. 63. 

I Birkhoff, loc. cit., p. 375, or Bdcher, loc. cit., p. 405. 
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conditions were such as to require a different solution of the equation between 
each two points of the set Oj , i = 1,2, • • ■ ,k. We can define a set of con- 
ditions which will do this as follows: Let wa be of the same form as Wij, 
but with different values for the a's, and let Wa be considered as applying 
to the solution to the right of the point a,, while wa applies to the solution 
to the left of that point. Then the auxiliary conditions are 

Wi{u) = Wa{u) + Wi^{u) + Wii{u) + • • • 

+ Wi,k-l{u) + Wi,k-l{u) + Wik{u) 
(i = 1,2, ■.■,n{k - D). 

When the general solution of i ( m ) = is substituted into these auxiliary 
conditions the constants are determined differently between each two of the 
points Oj. If that part of the condition Wi which applies to the interval 
(07, a-m) , i. e., Wii (ti) + Wim{u') , be denoted by Wu, then the determinant 
whose vanishing is the necessary and sufficient condition that the reduced 
system have solutions, other than the one which is identically zero through- 
out the interval ( a j^ 6 ) , becomes 



(4) Z) = 



W^ii(2/i) •■■ ^1,4-1(2/1; 

fF«(*-l), 1(2/1) ••• W^nCfc-l), 4-1(2/1) 

Fn(2/2) •••TFi.,_x(2/2) 

IF2i(2/-2) ••• WKo,i_l(2/2) 



W^«(i-i), 1(2/2) ••• W n(k-vi . k-\{yt) •■• 

i^ii(2/n) ••• W^,,_i(2/„) 
W^2l(2/») ••• W^2.*-l(2/n) 



Wn(k-V) , 1 ( 2/n ) • ■ • Wn(h-V, , i-1 ( Vn ) 



As in the case already considered, one easily finds that the necessary and 
sufficient condition that the complete system have one and only one solution 
is that the reduced system be incompatible. 

If it is assumed that the auxiliary conditions are linearly independent, the 
Green's function for the present system is defined for all values of x and s 
in ( a , 6 ) by requiring that for every value of s it shah satisfy as a function 
of X the following conditions: 
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I. It is continuous in the interval {a, b) , together with its first n — 2 
derivatives, except at the points at; 

II. Apart from these points, it satisfies L{u) = throughout the interval 
except at the point s, where its {n — l)th derivative has a finite jump of 
magnitude 1; 

III. It satisfies the auxiliary conditions Wi{u) = (i = 1, 2, •••, 
n{k - 1)). 

The proof that these three conditions completely define G{x, s) is parallel 
to that given for the simpler case. 

It is now possible to define adjoint auxiliary conditions, which, taken with 
the adjoint equation, form the adjoint system. The differential expressions 
L{u) and M (v) are connected by Lagrange's identity, 

vLiu) -uMiv) =-^Piu,v), 

where P is a non-singular bilinear form in the 2n variables u, u' , u" , • • • , 
y (n-i) . J, ^ ^' ^ ^" ^ . . . ^ ^(n-i) ^ for every x in the interval. If we integrate this 
over the interval ( o , 6 ) , we can for the purposes of integration divide it up 
into the sum of the integrals taken between each two successive points of 
discontinuity of 6? , i. e., the points aj , and thus admit the possibility of using 
solutions of the two equations with finite jumps at these points. This gives 
us 

(5) { [vL{u) -uM{v)]dx=^P, 

where P is now a bilinear form in the two sets of 2n ( fc — 1 ) variables each 

M«>(on, u^^{(h), u^Hat), u^'Has), ■■■, M«)(aLi), u^'Hau) , 

v^Hat), ««>(a2-), v^^Hat), v^^Hm), •••, v^^Hat-i) , v^'Hal) 

U = 0,1,2, •••,n-l). 

If Wi (u) ,i = 1,2, • • • ,2n{k — l),be any 2n{k — 1 ) linearly independent 
linear forms in the U's, then the U's can be expressed linearly in terms of 
them, and P becomes a linear form in the W's with coefRcieiits which are 
linear in the v's and which we may denote by Vt . So we can write 

P^Wi V2nCk-l) + W2 V2n(k-1)-1 + ■■■ + W2nik-1) Fj . 

In particular we may choose as the first n ( ^ — 1 ) of the linear forms in the 
u's the left-hand members of our n{k — 1 ) auxiliary conditions. Then the 
V's with indices 1 to n{k — 1) when set equal to zero give the adjoint auxiliary 
conditions. It is easily seen that different choices of the last n{k — 1 ) of the 
W's would lead merely to linear combinations of these adjoint conditions, 
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SO that the system thus defined is unique. From the symmetry in the work 
it follows that the adjoint relation is a reciprocal one. 

By definition G(x, s) as a function of x satisfied the reduced auxiliary 
conditions. We shall now prove that as a function of s it satisfies the adjoint 
auxiliary conditions. By substituting, in the integral and in P , G {x, s) 
for u and G{s , x) for v , we obtain 

''zV2nc.-iw(G)F,(G) = 0. 

The IF's here appearing are the last n{k — 1) . We shall make n{k — 1 ) 
different choices of these. We obtain in this way n{k — 1 ) different sets 
of V's, but as we have already remarked, any set can be expressed as linear 
combinations of the first set. If we do this we have n{k — 1 ) linear homo- 
geneous equations in the n{k — 1 ) variables, the first set of V's. The deter- 
minant of these equations depends on the choice of the W's, but it is easily 
shown that these can be chosen in such a way as to make the determinant 
different from zero, and it follows that the V's are all zero. That is, the 
Green's function as a function of its second argument satisfies the adjoint 
auxiliary conditions. Combining this with what we already know about 
G{x , s) as a function of its second argument, we see that (— 1)" G{x , s) 
is, as a function of s , the Green's function of the adjoint system. 

If yij (x) is defined as equal to yi{x) in the interval ( a.j , aj+i ) and as 
zero elsewhere in the interval {a, h) , then the n functions y and the k — 1 
intervals give in all n ( A; — 1 ) linearly independent solutions of i ( w ) , which 
it is convenient to call a fundamental system of discontinuous solutions with 
regard to the points a,. Such a system has properties analogous to the 
properties of a fundamental system of solutions in the ordinary case. Any 
solution of i (w) = having finite jumps at any of the points ai, but other- 
wise continuous in ( a , 6 ) , can be represented as a linear combination of these. 
We shall call any function which may have finite jumps at any of the points at , 
but satisfies L{u) = at all other points of ( a , 6 ) , a discontinuous solution 
(with regard to the points a,- ) . We have then the 

Theorem. If the reduced system has h linearly independent discontinuous 
solutions the adjoint system has also. 

Because of the reciprocal relation between the two systems it is sufficient 
to show that it has at least h . To prove this let y{x) represent any dis- 
continuous solution of i ( M ) = , and form a fundamental system of dis- 
continuous solutions, Zij (x) {i = 1 , 2 , • • • , n; j = 1 , 2 , • • • , k — 1) . Sub- 
stitution of these in the integral and in P leads to the following equations: 

n(A:-l) 
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X) W2n(k-l)-j {y)Vji{Zi) = , 



T. Winik-i)-j{y)Vj,k-l{Zi) = Q (i = l,2 ••-,«). 

These may be regarded &s n{k — 1 ) homogeneous linear equations to deter- 
mine the constants W2n(k-i)-j(y) , i = l, 2, 3, •••, n{k — 1) . Their 
determinant is just the determinant whose rank determines the number of 
Hnearly independent discontinuous solutions of the adjoint system, as will be 
seen by comparing it with (4), which is the corresponding determinant for the 
reduced system. We wish then to show that the rank of this determinant is 
at most r — h, which we do by showing that the system of equations has at 
least h independent solutions. To this end let us substitute the h linearly 
independent solutions of the reduced system in place of the solution y (x) 
in the W's. This gives h solutions of the system of equations, which can 
easily be shown to be linearly independent. So the theorem is proved. 

Of course these solutions which we have been dealing with may in particular 
cases be continuous over the whole interval (a, b) . We can regard our 
system (2) as a special case in which the first n auxiliary conditions are the 
conditions there given, and the rest of the set are the n{k — 2) conditions 
that the solution of the system together with its first n — 1 derivatives be 
continuous across each of the k — 2 points a, interior to the interval {a, b) . 
We have shown, then, how a system adjoint to the system (2) may be defined, 
and from now on the discussion will be limited to such a system and its adjoint. 

2. Characteristic numbers and solutions 
From the system (2) is formed the system 

(6) X(w)+Xm = 0, Wi{u)=0 (i = 1,2,3, ••■,n), 

in which X is a complex parameter with the whole finite plane for its domain. 
The adjoint system can be written 

(7) M{v)+'Kv = 0, Vi{i})=0 (i = 1,2,3, •••,n(fc - 1)). 

The solutions oi L{u) +Xm = are all functions of X , and so the determinant 
(3) is also. The values of X for which this determinant vanishes are the values 
for which the system above will have solutions other than the identically zero 
one. Such values are called characteristic numbers, and the corresponding 
solutions are called characteristic solutions. The characteristic numbers of the 
system and its adjoint are, by the theorem just proven, the same; and further- 
more the number of linearly independent discontinuous solutions corresponding 
to any one characteristic number is the same for both systems. 
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If u and V are solutions of (6) and (7) respectively, belonging to different 
characteristic numbers, then 






b 

u(x)v(x)dx = 0, 



as will readily be seen by substituting them in the integral (5). 

If we assume that each characteristic number has but one corresponding 
characteristic solution, and try to expand an arbitrary function in terms of 
these solutions, and if we assume further that the series is such that it can be 
multiplied through by a continuous function and integrated term by term, 
then we find that 



00 



(8) fix) = ^CiUiix), 

where 



4 



^» ^~ ^b 

I Ui{x)Vi 



f {x)Vi{x)dx 

*J a 

Ci = ,.„ 

{x)dx 



Now the Green's function of the system (6) is, of course, a function of X , 
and is defined for all values of X save the characteristic numbers. The funda- 
mental system 2/1,2/2,2/3, ■ ■ • , Vn, may be chosen analytic in X , whence 
G{x, s\ X) is also, except possibly at the characteristic numbers. If X = Xj 
is a characteristic number to which corresponds but one independent charac- 
teristic solution, and for which G {x , s; X) has a pole of the first order, then 
the residue at the pole is 

Ui{x)Vi{s) 



f 



Ui{x)vi{x)dx 



and the integral in the denominator is not zero. The proof given by Birk- 
hofl* applies here. 

Under these conditions the terms of the expansion (8) are 

(9) f Ri{x,s)fis)ds, 

in which Ri{x , s) is the residue corresponding to the ith characteristic num- 
ber, and the sum of the first m terms is 

(10) ~.jjy{x,s;\)f{s)dsdX, 

where F is a contour in the X-plane enclosing the first m characteristic numbers. 
* Birkhoff, loc. cit., p. 378. 
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In a previous paper* the author has proved the convergence of the integral 
(10) when suitable restrictions are placed on/(a;) and the auxiliary condi- 
tions. Of course in the more general cases the expansion (9) cannot be written 
in the form (8), that is, the residue cannot be expressed linearly in terms of 
the characteristic and adjoint functions alone. The new functions which 
must be introduced have been defined by means of integral equations by 
Goursat,t but apparently have never been defined for the two point case in 
terms of the differential system. Such a definition has been formulated by the 
author, and the accompanying theory forms the subject of another paper. 

* Wilder, loc. cit. 

t Goursat, Annales de la Faculte des Sciences de Toulouse, 
series 2, vol. 10 (1908), pp. 5-98. 



